Theory of the Metal-Paramagnetic Mott-Jahn-Teller Insulator Transition in A4C60 
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We study the unconventional insulating state in A4C60 with a variety of approaches, including den- 
sity functional calculations and dynamical mean-field theory. While the former predicts a metallic 
state, in disagreement with experiment, the latter yields a (paramagnetic) Mott-Jahn-Teller insula- 
tor. In that state, conduction between molecules is blocked by on-site Coulomb repulsion, magnetism 
is suppressed by intra- molecular Jahn- Teller effect, and important excitations (such as optical and 
spin gap) should be essentially intra- molecular. Experimental gaps of 0.5 eV and 0.1 eV respectively 
compare well with molecular ion values, in agreement with this picture. 
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The alkali fullerides A n Cgo, based on tightly bound 
Ceo molecules, are relatively simple molecular conduc- 
tors where strong correlations, orbital degeneracy and 
Jahn- Teller effect intervene at the same time, leading to 
remarkable physical properties jy. The three-fold de- 
generate ti u state of the isolated molecule leads to three 
bands hosting the valence electrons provided by the al- 
kali metals. Hence all compounds A n C6o with < n < 6 
should in principle be metallic if, as one expects, the 
bandwidth (of order 0.5 eV), is larger than the relatively 
small crystal field splittings. The n — 3 fullerides are 
in fact generally metals, and become superconducting 
with T c as high as 40 K Among the n = 4 com- 
pounds, however, Na4Cgo, a stable fee structure at high 
temperature ||], is the only one which is also a param- 
agnetic metal, whereas K4C60 and Kb^Ceo, stable bet 
structures down to low temperatures, are instead param- 
agnetic narrow-gap insulators. Understanding in detail 
this insulating state proves nontrivial, and constitutes 
the main scope of this work. Known experimental pa- 
rameters for K4C60 are a minimum band gap (probably 
indirect) between 0.05 and 0.2 eV and a direct, op- 

tical gap of 0.5-0.6 eV 0. The insulating state is non 
magnetic with a spin gap to the lowest triplet exciton of 
0.1-0.14 eV f6,||. Transition to a metal can be provoked 
by pressure |6|, which rationalizes why Na4Cgo should be 
metallic, Na being the smaller cation. 

One can invoke at least two possible scenarios for the 
insulating behavior of A4C60: a) a band insulator, due 
to a strong splitting of the t\ u bands, arising for example 
from the bet distortion, or else from a collective static 
Jahn- Teller (JT) distortion of the Ceo molecules, or al- 
ternatively b) a Mott-Jahn-Teller insulator ||, where 
the hopping between adjacent molecules is first of all 
suppressed by a strong Coulomb repulsion, an intra- 
molecular JT effect subsequently optimizing the state of 
the four localized electrons. We have studied first the 
molecular ion, and then the possible origin of the insulat- 



ing behavior of A4C60 by several theoretical approaches, 
tight-binding Hartree-Fock, density functional, and Dy- 
namical Mean Field Theory (DMFT), the latter a pow- 
erful method to treat the Mott transition ficfl . 

We start with the single CgQ~ molecular ion. Within the 
t\ u orbital, assuming rotational (icosahedral) symmetry, 
a general interaction among electrons can be written as 
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where n,i — Y~L c; „c, • is the electron number on each 
orbital (i = 1,2,3), n = n\ + 112 + TI3, and Ay — 
J2a c la c ja + H.c. At fixed n, the second term is re- 
sponsible for the multiplet splitting. For n = 2, 4, the 
lowest energy state is the 3 Ti g Hund's rule triplet, fol- 
lowed at energy 2U2 by a 1 H g singlet, and at 5U2 by a 
1 A g singlet. The exchange coupling U 2 is expected not 
to be significantly screened by the t\ u electrons, so that a 
reasonable estimate can be obtained by optical measure- 
ments on solid Ceo, which give U2 — 0.05eV |11 . Next, 



we consider the JT coupling to the eight H g vibrational 
modes. In the adiabatic limit, one can use a single-mode 
approximation ]l2] | 

Hjt = {{z +r ) + 

z (n\ + ri 2 — 2n 3 ) + rV3 (n\ — n 2 )| • (2) 

By using vibrational frequencies and couplings extracted 
from gas phase Cg photoemission [ p"3| , we obtain fro;, = 
0.117eV, and a dimensionless coupling g„ = 1.204, so that 
Ejt — h<->J*g1 = 0.169eV. Instead of attempting an exact 
solution of Hi n t + Hjt, we consider two opposite limits, 
antiadiabatic and adiabatic. Although neither of them 
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strictly applies, for hu* is comparable to all splittings 
and the coupling is of medium strength, the former limit 
will yield the right symmetries, while the latter will be 
quantitatively much more accurate. 

In the antiadiabatic limit, the JT term (Q) gives rise to 
a non-retarded electron-electron interaction, which can 
be absorbed into U2 1 — ► U2 — (3/4)£^jy, with a change 
of sign of U 2 from 0.05 to -0.076eV. The lowest en- 
ergy state is now the 1 A g singlet, followed at 0.23eV by 
the 1 H g singlet, and at 0.38eV by the 3 T lg triplet. The 
overall JT energy gain in this limit is very large, 0.84eV, 
about a factor three larger than the bare adiabatic JT 
energy (see below). This signals a true enhancement, 
due to the gain in zero point energy corresponding to 
the frequency collapse of the tangential vibron modes, 
first pointed out in Ref. H] , as a possible mechanism for 
explaining the high critical temperature of A3C60 com- 
pounds. 

In the adiabatic limit we diagonalize (0) plus (0) for 
n = 4, minimizing successively the eigenvalues with re- 
spect to z and r treated as classical variables. We find a 
lowest energy singlet at a classical distortion z = —1.987 
and r = 0, gaining Ejt — 0.293eV. Adding the zero 
point energy gain hu>* p4[ , we obtain a total gain of 
0.41eV, the zero-point enhancement still sizable. (A sim- 
ilar total gain of 0.42eV was obtained by uncorrelated 
eight-mode calculations [|l3|,^2).) The lowest triplet state 
has instead z = 1.0 and r = and lies above the ground 
state by E t — 0.108eV (spin gap). The lowest singlet, 
(with z = 1.0 and r = 0), is at E s = 0.208eV above the 
ground state. Finally, the optical gap, identified with the 
JT orbital splitting, is A = (3/2)zE JT ~ 0.504eV. 

Both limits therefore predict that the Cg^ ion should 
be a singlet. If hybridization between adjacent molecules 
were much smaller than all the molecular gaps involved, 
then a lattice of Cg^" molecules would indeed be a non 
magnetic insulator, moreover with an optical gap A and 
a spin gap E t remarkably close to the experimental ones. 
In that case, an electronic structure calculation and total 
energy minimization for the A4C60 lattice should yield 
a narrow t\ u band split by an insulating gap, in turn 
supported by a static (uniform or staggered) collective 
JT distortion of all molecules (scenario (a)). 

A tight-binding Hartree-Fock (HF) approximation 
yields precisely that. The Hamiltonian which we use is 

H= ^Ha^ + Hint + HjT, (3) 

where the hopping amplitudes tfj are evaluated along 
the lines outlined in Ref. jlj| , and the vibronic terms are 
treated in the adiabatic approximation, as in Eq.(§). We 
look for a HF state with a non-zero uniform average of 
m = n-2 ^ ri3 (collective JT state) and find a stable band 
insulator with direct gap and spin gap of ~ 1.48eV and 
an indirect gap of ~ 1.03eV, all of which are much larger 



than the experimental values. The calculated HF band- 
structure is shown in Fig. 1 (upper panel) . This indicates 
that this method, known to overestimate insulating ten- 
dencies, is unreliable for this problem, and more realistic 
first principles calculations are called for in this case. 

We have carried out a series of such calculations for 
K4C60, starting with bet configurations compatible with 
X-ray data [|l6| , and using state-of the-art pseudopoten- 
tial plane- wave density functional techniques [jl7|], with 
E cu t = 35 Ry (55 Ry for refinements), and careful /c-point 
summations. Confirming previous results [Q, we found 
first of all that intermolecular electron hopping is not 
small, yielding metallic ti u bands of width W ~ 0.6eV. 
We searched next for a spontaneous collective JT dis- 
tortion by relaxing atomic positions based on Hellmann- 
Feynman forces. In order to check that the calculation 
could in principle yield such a delicate JT distortion, and 
also reproduce the molecular limit, we performed test cal- 
culations at an artificially enhanced intermolecular spac- 
ing of 13. 3A. In that case a distortion appeared sponta- 
neosly, leading to a distribution of carbon distances from 
the Ceo center between 3.511 and 3.553A. This distor- 
tion magnitude AR = 0.042A, although small, is very 
close to that reported for [PPN(+)] 2 C 60 (2-) salts |§, 
namely 0.043 A, and that is very gratifyng since JT dis- 
tortions of Cg 2 and C 60 7 should be essentially the same. 
The single-particle gap was only about 0.1 eV, instead of 
the expected 0.5 eV, a standard density functional short- 
coming of no consequence for this case. When carried out 
for the true structure (bet, intermolecular spacing 9.97A) 
however, an extensive search for a JT distortion with 
one molecule/cell yielded no result, and the system re- 
mained undistorted and metallic, as shown in Fig. 1 
(lower panel). 

Notably, X-rays fail to find a static JT distortion jl6) , 
however within a resolution of precisely 0.04A, which 
is inconclusive. Selected trial calculations with two 
molecules/cell also failed to yield a doubling, as would 
be caused, e.g., by charge-density- waves jl8| or by a stag- 
gered collective JT state. We temporarily conclude that, 
since accurate density functional calculations cannot ac- 
count for its insulating behavior, K4C60 is probably not 
a band insulator. We will make this statement stronger 
based on additional reasoning at the end of this paper. 

Electron correlations can be strong in these com- 
pounds. A realistic estimate p]Jll|, of the intra-molecular 
Hubbard U of Eq.(|l|) is ~ 1.0-1.6eV, which is larger than 
the full bandwidth. The failure of density functional cal- 
culations does not appear surprising in this light. To get 
a description of the insulating state, and of the insulator- 
metal transition, we resort to dynamical mean-field the- 
ory (DMFT) Jl0| , which is exact in the limit of infinite 
coordination lattices, and has proved quantitatively suc- 
cessful in describing the Mott transition. 

Within DMFT, the Hamiltonian ([|) is mapped onto 
a three-fold degenerate impurity Anderson model, where 
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the parameters describing the conduction electrons, as 
well as their coupling to the impurity, are determined self- 
consistently |lC| , Band structure enters the calculations 
via the density of states (DOS) in the self-consistency 
equations. We approximate the realistic DOS with a 
semicircular DOS of same bandwidth W, appropriate for 
a tight-binding model on an infinite-coordination Bcthe 
lattice. This approximation does not change the quali- 
tative behavior of the model and is not expected to sig- 
nificantly influence the value of the gap. The Anderson 
model is solved by means of exact diagonalization with 
a finite number (n s ) of conduction electron degrees of 
freedom |^of , checking convergence as a function of n s 
by finite-size scaling. We should look for either orbital 
symmetry broken and unbroken phases, as the capability 
to describe a true Mott insulator without any symme- 
try breaking is a unique feature of DMFT, which is not 
shared by HF or density functional-based calculations, 
always implying symmetry breaking at metal-insulator 
transitions. For simplicity, we have restricted our search 
for insulating solutions either without orbital symmetry 
breaking or with uniform JT ordering. In order to locate 
the metal-insulator transition (MIT), we calculate the 
quasiparticle residue Z (plotted in Fig. 2 as a function of 
U/W), whose vanishing identifies the critical U /W above 
which a paramagnetic insulator is stable, and the metal- 
lic one is not. Again, we limit our analysis to limiting 
cases: (i) No JT coupling; (ii) Antiadiabatic JT effect; 
(Hi) Adiabatic JT effect. Case (i) - Orbital degeneracy 
without JT coupling, is well studied for C/ 2 =0 Jl|,^l|. 
In the presence of a finite U 2 /W = 0.08 we find a MIT 
at U/W ~ 1.414 for symmetric solutions, much reduced 
with respect to U/W ~ 1.98 with U 2 = 0. Case (ii) - The 
JT coupling renormalizes the dipolar integral, leading to 
an effective U 2 /W = -0.127. The critical MIT U/W is 
shifted to a much lower value U/W ~ 0.707. Case (Hi) 
- Solution of DMFT equations is more involved, since 
the self-consistency must be required only after averag- 
ing over the classical vibrational (five-dimensional) vari- 
able Q 0,E3|. In the broken-symmetry case the aver- 
age simplifies, at T = 0, as the probability distribution 
of Q becomes a single 5-function. We obtain a critical 
U/W ~ 0.9 — 1.0 (the uncertainty due to convergence 
difficulty). The orbitally symmetric insulating solution, 
which becomes stable at U/W ~ 1.237, is of particular in- 
terest since it describes a molecular insulator where each 
molecule is distorted with equal probability in all possible 
directions and independently from any other molecule. 
This state has therefore a very large entropy which could 
be reduced by including quantum fluctuations, for in- 
stance in the form of tunnelling between the equivalent 
local distortions, a way of describing a dynamical JT ef- 
fect, not included at the adiabatic level. It should be 
noted that our adiabatic limit does not now include the 
vibron zero-point energy gain. Because it is present in 



the symmetric case only, it could lower the true critical 
U/W value of this phase, which could in reality prevail 
over the broken symmetry case. Moreover, temperature 
would also favor the symmetric state, where entropy is 
higher. 

We can now compare critical U /W's with the true ones. 
As before, we expect the adiabatic values to be quantita- 
tively more accurate. The calculated critical values are 
in all cases substantially smaller than the actual U/W 
value, i.e. 1.7-2.7. We conclude that the insulator is best 
explained as a Mott- Jahn- Teller state, where orbital de- 
generacy has becomes split, giving rise to an essentially 
intra-molecular Cqq Jahn- Teller state whose calculated 
gaps (optical, spin) agree very closely with experimental 
data. In the adiabatic limit the orbital symmetry is bro- 
ken, but that is likely to change when zero-point energy 
and temperature are included. 

The close agreement between the adiabatic gaps of Cqq 
and experimental optical and spin gaps in A4C60 is per- 
haps the strongest piece of evidence in favor of a Mott 
state. Suppose one could even find, by some other band 
calculation, such as GW H|, a stable static collective 
JT state as in scenario (a). By necessity, the collective 
JT distortion magnitude would have to be substantially 
smaller than that of the isolated molecular ion, since 
electrons leaving the molecule very frequently to hop on 
other molecules weaken the on-site JT effect. But, in 
that case, it should not be possible to observe optical 
and spin gaps of exactly the right molecular magnitude, 
as one does. They would be much smaller, correspond- 
ing to the derealization, or spillout, of the band Wan- 
nier function. In the Mott state, the electron spillout to 
neighboring molecules is reduced to order t/U, which is 
very small. Hence this is the only state that can explain 
why JT electronic gaps are essentially intra-molecular. 
When the Mott insulator state is destroyed, for exam- 
ple by pressure-induced increase of W and decrease of 
U, this intra-molecular physics is expected to disappear 
rather suddenly (giving way to a density functional-like 
metal), instead of gradually as in a band MIT transition. 

Our conclusion is that the bet A4C60 are not conven- 
tional band insulators, and that both strong correlations 
and JT effects are crucial for their understanding. That, 
incidentally, suggests that also the n=3 superconducting 
fullerides, whose bandwidth is quite similar, are most 
likeky close to a strongly correlated state |n]]. In that 
case the metal competes in our picture with a Mott-Jahn- 
Teller spin 1/2 antiferromagnetic insulator, with the im- 
plication that strong correlations should probably not be 
ignored when discussing superconductivity. 
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FIG. 1. The tin bands for K4C60 according to 

Hartree-Fock (upper panel) and LDA (lower panel). The HF 
state sustains a collective JT distortion and is wide gapped, 
the LDA state is undistorted, and metallic. Both approxi- 
mations disagree with experiment, indicating a narrow gap 
semiconductor. 
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FIG. 2. Quasiparticle residue Z as a function of U/W 
for the symmetric solution of the purely electronic model 
(solid line), the symmetric (dotted line) and broken symme- 
try (dot-dashed line + heavy dots) solutions in the adiabatic 
limit, and the symmetric solution in the antiadiabatic limit 
(dashed line). The adiabatic broken symmetry value should 
be the best numerical approximation. Realistic values of U /W 
clearly indicate a Mott-Jahn- Teller insulator. 
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